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Abstract 

Evolution of structure functions is studied in a chiral potential model which incor- 
porates the pion as the chiral symmetry restoring field. Evolution equations for 
the quark and pion densities are derived in the lowest order at large Q^. The split- 
ting function for quark emission from point-like pion is flavor-dependent. This is 
shown to lead to nontrivial evolution of the nonsinglet moments in the next order 
which is consistent with the observed departure from the Gottfried sum rule. 
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1 Introduction. 

The observed departure from the Gottfried sum rule is difficult to explain 
within the framework of perturbative QCD [1]. The effect indicates a flavor- 
asymmetry in the sea quark distribution inside a nucleon which cannot be gen- 
erated by the Altarelli-Parisi evolution of the structure functions. The simplest 
alternative is to bring in pions into the picture [2,3,4] and replace gluon emission 
from a quark by pion emission [5]. Just having pions inside the nucleons cannot, 
of coiuse, explain the observed effect because the direct contribution of the pion 
density drops out from the equation of evolution of the nonsinglet moment which 
is the object of interest in the context of the Gottfried sum. A rather complete 
treatment of the problem has been given by Ball and Forte [6] who have com- 
puted the non-perturbative cross sections for the emission of pseudoscalar bound 
states from quarks and obtained a satisfactory explanation of the observed phe- 
nomenon. Their work shows that the effect is essentially nonperturbative. Also, 
there is good theoretical ground to expect such evolution, even if the analysis of 
the experimental data is clouded with uncertainties. 

In this paper we are interested in seeing whether such an evolution can result 
in the phenomenological models of the nucleons which have been devised to study 
the static properties of the nucleon bound state. The question is nontrivial since 
the connection between these models and the QCD-evolved structure functions 
is not at all obvious. On the other hand ,the models of the nucleon bound state 
are supposed to be phenomenological manifestations of the nonperturbative QGD 
at the hadronic scale[16]. It is therefore natural to try to see to what extent 
the models serve as effective theories and whether they can be put to any use 
other than calculation of static properties of the nucleon. As an example, we have 
chosen in this paper a chiral potential model incorporating point-like pions and 
have calculated the relevant cross sections and the splitting functions using the 
quark-pion coupling of that model. The choice of this particular model is not due 
to its having any greater intrinsic merit than other existing models, but due to 
the simplcity of its quark-pion hamiltonian. 

We report here only the low order perturbation results at large and show 
that a nontrivial evolution of the nonsinglet distribution can indeed result from 
this model. 

2 The model. 

We consider, with some modification, a phenomenological independent-quark 
model of baryons studied by Jena and Panda[6,7] (sec also [8-15]) who could ob- 
tain the static electromagnetic properties of the nucleon, the axial vector coupling 
constant in beta decay and the pion-nucleon coupling constant in resonable agree- 
ment with the experiments. The quarks in this model move independently in a 
confining potential Vq{r) which is an equal mixture of a vector and a scalar part: 

V,{f^ = lil + r)V{r) (2.1) 
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with 

V(r} = a^r + Vo, a> o 

The non-invariance of the quark mass term and the scalar part of Vq{r) under 
SU(2) chiral transformation is sought to be adjusted by introducing a pion field 
interacting with the quarks through a linearised term in the lagrangian density 

£]{x) = f G(f)V^,(x)7V.0(x)V;,(x) (2.2) 

G{r-) = m, + ^ (2.3) 

TTiq is the quark mass and /tt is the pion decay constant. The quark wavefunction 
of the baryonic core is determined fromV^(r);in the ground state configuration 

\ Xq dr / 

(t>q{r) turns out to be expressible in terms of the Airy function. 

The static properties of the baryons calculated fromVq(r) alone (with a ~ 
0.343 Gev,Vo ^ —0.506 Gev)are alredy close to the observed values so that, on a 
phenomenological level, the quark-pion coupling is expected to be small making 
perturbative calculations meaningful. On the other hand, this is an example of 
a potential model which allows one to calculate pion emission and absorption by 
quarks and their effect on the structure functions. We show this by considering, as 
the custom is, a quark beam which is being probed by a virtual photon, obtaining 
the splitting functions from the lowest-order cross sections. 

The relevant cross sections are (a) (j{j*q q't^), (b) cr{^*q —>■ 7rg')and (c) 
C7(7*7r qq). 

The processes (a) and (b) contain emission of quarks from quarks and pi- 
ous, and lead to the splitting functionsPgg andP^Tr-The process(c) gives the splitting 
functionPTrg- 

The external potential G having a long range part in the original model of 
Panda and Jena necessitates the introduction of an infrared cutoff parameter A in 
the three momentum of the static gluon. Physically, this is fixed by the size of the 
baryons. With this modification, the effective qqir vertex factor for perturbation 
calculations in this model becomes 



rqq7r(P,y) = 



1 



4a2(|p-p'| - A) 



7rA|p — 



7' (2.5) 



where p and p' are the quark four-momenta going into and out of the vertex. 
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3 Virtual photon cross sections in the large limit. 

The evolution of quark and pion densities in this model can be obtained by 
first calculating the virtual photon cross sections listed in the previous section. The 
lowest order perturbation calculation using the vertex factor(2.5) is straightforward 
in the small transverse momentum and large limit neglecting the quark and 
pion masses. 

In the following p, p',k,q denote the (center-of-mass)four-momenta of the 
quarks q, q' , the pion and the photon.For the process(a),we define the invariant 
variables Sa = {q + p)^, ta = {p' — q)^, Ua = {p — p')^- Then in the leading order, 
the cross section cr{^*q —>■ q'n) is 



da 



ttA 



(3.1) 



Defining Sh — {q + p)'^,tb — {q — and ui) — {p' — q)'^, the cross section for 
the process 7*9 —> 7rq' is 



da 

dnZ 



27ra 



ttA 



-^6 



(3.2) 



Defining Sc= {q-\- k) , tc = (q — p') and Uc = (p' — k) , we find, to leading 
order in Q^, the cross section for the process 7*7r — > q'q 



da 

dn' 



27ra e^, 



4a2 
ttA 



1 
2 



tr.Ur 



- -1 



(3.3) 



To facilitate the extraction of the splitting functions, we express these cross sections 
in terms of the kinematic variables z and pt for each process which are related to 
the momentum variables by 



z ~ 



p% ~ 



-St 



the relations being true for small t and large Q^. Under these approximations, 

'l + z^ 



(3.4) 
(3.5) 



da{pi*q q'n) _ 2 

^ Pt 



dp 



1-z 



da{-f*q ^ nq') 2 "tt^ 
— = e^a-o-2-22; 



Pt 



da{-f*TT q'q) _ ^2 



dp: 



— C^/CTq 2 



Pt2 



{z^ + {l-z)')+[^\z{l-z) 



(3.6) 
(3.7) 

(3.8) 
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where Uo = Att a/s and 



1 



14 



4a^ 
ttA 



(3.9) 



which is an effective couphng constant of the model in the hmit of large Q^. 

4. The evolution equations and splitting functions. 

Integrating (3.6)- (3.8) with respect to (with an infra-red cutoff), we can 
now easily get the evolution of the structure function F2(a;, Q^). The quark den- 
sities are then found to satisfy the following equations. 



= a^ f ^ [u{y)Puu{z) + d{y)Pud{z) + 7r°{y)Pu^o{z)+7r+{y)P^^+{z)] 

J X y 

(4.1) 

= a^ [ ^[u{y)Puu{z)+d{y)P^aiz)+'^°{y)Pu.^{z)+7v-{y)P^^^{z)] 

J X y 

(4.2) 

and two similar equations for d(x, Q'^) and d{x, Q^). Pab stands for the splitting 
function for 6 — > a and z = x/y. The pion densities are found to satisfy 



du{x) 
dlogQ'^ 

du{x) 
dlogQ"^ 



d7r°{x) 
dlogQ"^ 



OL^ 



dix^ix) 
dlogQ"^ 

d7r~{x) 
dlogQ^ 



a J '^[u{y)P^^^{z) + d{y)P^^M] (4.3) 

Jo y 
dy 



[d{y)P^-d{z) + u{y)P^-^{z) 



(4.4) 



[u{y)P^o^{z) + d{y)P^oa{z) + u{y)P^o^{z) + d{y)P^od{z)] 

(4.5) 



The splitting functions are given by 



z 1 — z 



P,. = 2 



(.2 + (l-.)2)+(^)%(l-.) 



(4.6) 



The splitting functions obviously satisfy the requirement of charge conjugation 
invariance : Pqq — Pqq, PqT^+ = Pqjr- , Pq-n° — Pqn° etc.As given in (4.6), they 
do not satisfy the fermion number conservation, but this could be ensured by the 
usual regularization procedure. 

5. Evolution of the nonsinglet distribution. 

The first thing one notes about the results given above is that even in the low- 
est order the evolution of individual quark densities is flavor-dependent. However, 
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this alone can not explain the observed evolution of the Gottfried sum. Firstly, the 
first order calculation of the cross sections will always lead to Pqq = Pqq = O.The 
contribution of Pg/^r drops out from the evolution of the nonsinglet distribution, 
and hence the dependence of P^tt does not matter at all in this context. 

Secondly, the drastic large approximation used to obtain the results of 
the previous section will not lead to any evolution of the anomalous dimensions 
obtained by taking moments of the splitting functions. Hence, the Gottfried sum 
remains unaffected as it is controlled by the anomalous dimension for the first 
moment. 

Therefore, in this model it is essential to consider both higher order pertur- 
bation terms and low corrections. The splitting function Pqq is generated by 
the second order process and is seen to be of the form 



and gives a contribution to the anomalous dimension of the first moment 



This then points to a detectable evolution of the nonsinglet first moment and, 
hence, of the Gottfried sum.The anomalous dimension (5.2), obtained by retaining 
the large approximation, is not dependent on and, therefore, cannot conform 
to the observed departure from the Gottfried sum rule. The dependence appears 
when one considers the low corrections. 

6. Conclusion. 

The simple chiral potential model considered here leads to a nontrivial evo- 
lution of the nonsinglet moments even in the large approximation. This is a 
nonperturbative phenomenon in QCD. The phcnomenological model arrives at it 
by allowing for emission of an isovector meson by a quark. The detailed results 
for the evolution of the Gottfried sum with low Q"^ corrections will be reported 
elsewhere. 

One of us (K.R.-M.) acknowledges financial support from the Council of Sci- 
entific and Industrial Research during this work. 




TT 




(—\ -2\{z{l-z)+zHogz) (5.1) 




(5.2) 
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